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PQ . ABSTRACT 

^ \ Blazars constitute the most interesting and enigmatic class of extragalactic 

I \ 7-ray sources dominated by non-thermal emission. In this Letter, we show how 

;_i \ the WISE infrared data make possible to identify a distinct region of the [3.4]- 

^ \ [4.6]-[12] /im color-color diagram where the sources dominated by the the thermal 

radiation are separated from those dominated by non-thermal emission, in par- 
ticular the blazar population. This infrared non-thermal region delineated as the 
WISE blazar Strip (WBS), will constitute a new powerful diagnostic tool when 
O \ the full WISE survey data is released. The WBS can be used to extract new 

Q I blazar candidates, to identify those of uncertain type and also to search for the 

^ I counterparts of unidentified 7-ray sources. We show one example of the value 

O ■ of the use of the WBS identifying the TeV source VER J0648+152, recently 

discovered by VERITAS. 

^^ I Subject headings: galaxies: active - galaxies: BL Lacertae objects - radiation 

^ ' mechanisms: non-thermal 






1. Introduction 

Blazars are one of the most enigmatic and rare class of active galactic nuclei (AGNs). 
Their continuum emission is dominated by non-thermal radiation from radio to 7-ray ener- 
gies, making them the most frequently detected class of extragalactic sources at GeV-TeV 
energies. Their observational features also include high and variable polarization, super- 
luminal motions, very high observed luminosities coupled with a fiat radio spectrum that 
steepens in the IR-optical bands and a rapid variability from the radio to X-ray bands with 
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weak or absent emission lines. In 1978, Blandford and Rees suggested that radiation of 
blazars can be described as arising from a relativistic jet closely aligned to the line of sight 



( iBlandford fc Rees 19781 ). In short, blazars constitute the purest non-thermal kind of AGNs. 



The spectral energy distributions (SEDs) of blazars include two main components: a 
low-energy component with power peaking in the range from the IR to the X-ray band, and a 
substantial high-energy component often dominated by 7-rays. Blazars come in two flavors: 
the BL Lac objects (BZBs) and Flat Spectrum Radio Quasars (BZQs), with the latter 
having strong emission lines, generally higher radio to optical polarization higher redshift 
and more prominent 7-ray bump of their SEDs, than the former population. The BL Lac 
population is in turn divided in two subclasses: the "low-frequency peaked BL Lacs" (LBLs) 
in which the peak of the flrst component falls in the IR-optical range, and the "high-frequency 
peaked BL Lacs" (HBLs) w hen it falls in the range between the UV and the X-ray bands 



( jPadovani fc Giommi 19951). Recently, HBLs detecte d at TeV energies have been reclassifled 



as TBLs (JMassaro et al. 20081 : iMassaro et al. 201 laf ). 



In this Letter we present a color-color diagram using infrared (IR) magnitudes that 
allows us to distinguish between extragalactic sources dominated by non-thermal emission, 
like blazars, and other classes of galaxies and/or AGNs. We construct th e diagram using the 



archival data of the recent Wide Infrared Survey Explorer ( WISE) facility (IWright et al. 2010[ ) 
The WISE mission mapped the sky at 3.4, 4.6, 12, and 22 ^m in 2010 with an angular res- 
olution of 6".l, 6". 4, 6". 5 & 12". in the four bands, achieving 5a point source sensitivities 
of 0.08, 0.11, 1 and 6 mJy in unconfused regions on the ecliptic, respectively. Regarding the 
WISE astrometry we note that the absolute (radial) differences between WISE source-peaks 
and "true" astrometric positions anywhere on the sk y are no larger th an ~ 0".50, 0".26, 
0".26, and 1".4 for the four WISE bands, respectively JCutri et al. 201l[ l 0. 



We show how in the [3.4]-[4.6]-[12]/xm color-color diagram the blazar population, which 
is dominated by non-thermal emission in the infrared, covers a distinct region well separated 
from the locus of other extragalactic sources dominated by the IR thermal radiation. We 
suggest how this infrared non-thermal region (hereinafter the WISE blazar Strip, WBS) 
could be used as a diagnostic tool to identify the blazar population over the entire sky when 
the flnal data release of the WISE survey will be available. Finally, we discuss the possible 
implications of our results on efforts to identify 7-ray sources and blazar of uncertain type. 
As an example, we applied the WBS as a new diagnostic tool based on the WISE blazar 
Strip to look for the counterpa rt of the recently unid entifled BL Lac candidate discovered by 



VERITAS: VER J0648+152 (lOng fc Paneque 2010f ). We note that the WISE 22/im band 



^ http://wise2.ipac.caltech.edu/docs/release/preliin/expsup/sec2JBg.htinl 
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also offers a similarly useful color diagnostic, but because of its larger beam it is less useful 
in this investigation. 



Blazar WISE detections 



We considered all the blazars present in the ROMA-BZCAT § jMassaro et al. 2009 



Massaro et al. 20101 ). adopting the same nomenclature, namely: BL Lac objects (BZBs), 
fiat spectrum radio quasars (BZQs), and blazar of uncertain type (BZUs). We note that 
the class BZU which contains sources showing occasional presence/absence of broad spectral 
lines or transition objects between a radio galaxy and a BL Lac, as well as sources with 
limited spectral data does not allow a precise classification. 

The total number of blazars in the ROMA-BZCAT that fall in the area surveyed by 
WISE during the first year (corresponding to 57% of the whole sky) is 1487 (due to the 
inhomogeneity of the sky coverage of the ROMA-BZCAT). To search for the positional 
coincidences of blazars in the observed WISE sky we considered two different regions. The 
first of radius 2". 4 corresponding to the c ombination of the err or of 1" assumed for the radio 



position reported in the ROMA-BZCAT flMassaro et al. 20091 ) and that of the fourth WISE 



band at 22/im (i.e., 1".4) (see Section[T]and reference therein); The second region has a radius 
of 6". 5, equal to the angular resolution of the third WISEhand. In our analysis, we excluded 
multiple crossmatches. We also note that the cross-match between the ROMA-BZCAT and 
the WISE catalog has been performed considering only the sources present in the first year 
WISE catalog detected with a minimum signal-to-noise ratio (SNR) higher than 7 in at least 
one band. 

The number of positional coincidences within the first region of 2". 4 is 1365 correspond- 
ing to 92% of the blazars in the WISE first year survey, detected with a chance probability 
of 2.5%, evaluated applying the same method described in Maselli et al. (2011). We did not 
find any multiple matches adopting the area of 2.4". The number of blazars associated with 
WISE sources increases to 1446 with a chance probability of 11% considering the second 
region of 6". 5, where we found 8 multiple matches. Thus, 97% of the blazars appear to have 
a counterpart in the WISE catalog. To be more conservative in our analysis, we consider, 
for the rest of the discussion, only those blazars with WISE sources associated within the 
2". 4, unless stated otherwise. 



' http://www.asdc.asi.it/bzcat/ 
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3. The WISE hlazar Strip 



In the WISE Prehminary Source Catalog (WPSC) the [3.4]-[4.6]-[12] /iin color-color 
diagram, drawn fr om high and low galactic latitude regions, shows the location of different 



classes of objects (jCutri et al. 201ll ) 



In order to put blazar populations on the same plot, we selected 14 random regions of 
4 deg^ each, not overlapping, for a total 56 deg ^ at high Galactic latitude, within the 116 



deg^ considered in the WSPC (ICutri et al. 201 ll ). Then, we collected all the 453420 sources 



detected by WISE in its 1st year catalog (hereinafter WISE thermal sources, because they 
are dominated by thermal emission in the infrared), within the selected 56 deg^ region, 
having SNR>7 at le ast one band, a co nservative level for the WPSC release to emphasize 



Catalog reliability (ICutri et al. 20111 ). Being far from the Galactic plane, the majority of 



these sources have extragalactic origin and only few stars lie in the selected region. 

We built [3.4]- [4.6]- [12] //m color-color diagram using the magnitudes reported in the 
WISE Catalogs for all the sources present in the 56 deg^ and for all the blazars detected by 
WISE . We note that the relative errors for both the infrared colors are less than 0.1 for 
95% of the entire blazar sample but less to 0.05 for ~ 85%. In Figured! the two subclasses 
of blazars and those of uncertain type are shown, overlaid to the isodensity contours for 
all the sources in the 56 deg^ of the WISE selected area. The four levels of the isodensity 
contours correspond to 50, 100, 500, 2000 sources, in the 56 deg^ area, respectively. Finally, 
in Figure [H we also annotated the location of different classes of objects. 

The main striking result is that the blazars lie in a distinct region (i.e., the WBS) of the 
[3.4]-[4.6]-[12] /im color-color plot with respect to the WISE thermal sources, in particular 
where the density of the other extragalactic sources is decreasing that of the blazar population 
increases (see below for more details). Moreover the two main classes of BZQs and BZBs are 
themselves separated in two distinct regions of the WBS, where BZQs are generally redder 
and lie closer to the normal quasars (QSRs) and Seyfert galaxies. 

In Figure [21 we consider only the BZBs distinguishing between their subclasses of HBLs 
and LBLs. We calculated the ratio between the X-ray and the radio flux $xr for those 
source with these fluxes available in the ROMA-BZCAT to distinguish between HBLs and 
LBLs in the above sample. W e adopted a very c onservative classification considering as 



HBLs the BZBs with ^xr > 1 (IMaselh et al. 20101 ) and LBLs for all the others 



http://wise2.ipac.caltech.edu/docs/release/preliin/preview.htinl 



*A11 WISE magnitudes are in the Vega system. 
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Fig. 1. — The [3.4]-[4.6]-[12] /iin color-color diagram of WISE thermaA sources and blazars. 
We report the 1365 blazars associated with a WISE source within a region of radius 2". 4. 
The two blazar classes of BZBs (blue filled circles) and BZQs (red filled circles) are shown 
together with the blazars of uncertain type (BZUs, green filled circles). The background 
grey dots correspond to the 453420 WISE thermal sources detected in a region of 56 deg^ at 
high Galactic latitude. The isodensity curves for the WISE thermal sources, corresponding 
to 50, 100, 500, 2000 sources in the 56 deg^ area are reported (see Section [3]). The location 
of different classes of objects is also shown, where QSRs, ULRIGs and LIRGs indicate the 
quasars, the ultraluminous infrared galaxies and the luminous infrared galaxies, respectively. 
The WBS is highlighted within the two black dashed lines. 
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All the BZBs without X-ray data available are omitted from Figure |2] because it is not 
possible to classify them. The B ZUs have been excluded from our investig ation, because of 
the uncertainties in their nature (JMassaro et al. 20091 : iMassaro et al. 2010[ ). 



As many authors have pointed out the infrared spectrum of BL Lac objects, being dom- 
inat ed by non-thermal emission, can be descri bed as a power-law with the spectral index a 



[e.g. Ilmpey et al. 19821 : iBersanelli et al. 19921 ) {^Sy oc z/~°'). Then, diffe rent values of a cor - 
respond to different colors in the [3.4]- [4.6]- [12] /xm color-color diagram ( IWright et al. 2010l ). 
In Figure m we show the line corresponding to the power-law spectrum for different a values, 
ranging between 2 and -2. It is evident that the sources in the WBS are in agreement with 
the path described by the power-law model, confirming the non-thermal origin of their in- 
frared emission. There is also a marginal evidence that HBLs tend to lie below the 
power-law line model in Figure [2] while LBLs lie above the same. This deviation 
could be due to a contribution of the host galaxy or to a mild curvature of their 
infrared spectra in the WISE bands. Su ch interesting effect w ill be investigated 
in more details in a forthcoming paper ( IMassaro et al. 201 Ibl ) 



BL Lac objects 



HBLs 
LBLs 
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Fig. 2. — Same of Figure [Tj only for the BL Lac objects. The two subclasses of HBLs 
(black circles) and LBLs (red circles) are highlighted. We also report the line PWL (black) 
correspondent to an IR non-thermal power-law spectrum of spectral index a. We marked 
the IR colors of different a values in blue along the black dashed lines, (see Section [3] for 
more details). 
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Finally, we considered the list of HBLs detected at TeV energies (TBLs) as reported in 
the TeVcat |f|, updated on 2011 August, with ^xr > 1- Only 13 TBLs are present in the 
sky mapped by WISE during the first year, and all of them have a counterpart in the WISE 
archive within a region of radius 6". 5 while this number decreases to 11 if considering the 
smaller region of radius 2". 4. We marked the TBLs in Figure [3] in order to investigate if 
they have a different behavior from the TeV und etected HBLs (UBLs) as found in the X-ray 



analysis of a subsample of the ROMA-BZCAT ( jMassaro et al. 2011al ). We found that the 



TBLs are more concentrated near the center of the WBS (see Figure E]). 

To highlight the separation between the WBS occupied by the BZBs in the [3.4]- 
[4.6]-[12] //m color-color diagram and the region spanned by the overall distribution of 
WISE extragalactic sources, we have evaluated the two-dimensional densities of both pop- 
ulati ons of source using the Kernel Density Estimation (KDE) technique (Figure H]) (see 



e.g., iD'Abrusco et al. 20091 : iLaurino fc D'Abrusco 201ll . and reference therein). The KDE 
method provides an effective way of estimating the probability function of a multivariate 
variable and do not require any assumption about the shape of the "parent" distributions. 
In Figure HI the isodensity contours drawn from the KDE density probabilities and associ- 
ated with different levels of density are plotted for the two classes of sources (i.e., BZBs and 
WISE thermal sources). The separation between the region of higher density of the BZBs 
distribution and that of the general WISE source distribution is evident. 

In particular, the line A shown in Figure H] has been obtained by a linear regression of 
the BZB density population in the [3.4]- [4.6]- [12] /im color-color diagram, using the 95% of 
BZBs detected by WISE in 2". 4. It represents the main WBS axis. The second WBS axis 
(i.e., line B in Figure Hj) displays the direction of maximum gradient of the density ratio for 
the distributions of the two populations (i.e., BZBs and the VF75'£' thermal sources). The 
strong density gradient, visible along the line B, suggests a fairly good separation between 
the two distributions of BZBs and WISE thermal sources. The ratio of the BZBs density to 
the WISE theTmal sources density increases toward its maximum (i.e., point C in Figure H]) 
and along the curve A, while the density of the WISE thermal sources decreases. 

The separation of the WBS from the distribution of WISE thermal sources is also 
evident in the normalized density histograms of the distances |D^| from the line A for both 
populations (see inset of FigureH]). The peak of the BZB density histogram (inset of FigureH]) 
suggest a strong clustering of this source population, while that of the WISE ther-mal sources, 
with a completely different shape, peaks at a significantly large distance from the axis of the 
WBS . 



^ http://tevcat.uchicago.edu/ 



BL Lac objects 
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Fig. 3. — Same of Figure [H only for the BL Lac objects (blue circles). The TBLs are 
highlighted (green diamonds). The purple open squares indicated all the sources associated 
to the VER J0648+152 unidentified object within the error box of the Fermi position. The 
purple filled square points to the WISE source associated with the X-ray position of XMM- 
Newton, also corresponding to the second closest source within 0.03 deg. 
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Fig. 4. — The [3.4]-[4.6]-[12] fim color-color diagram for the BZB population and the WISE 
thermal sources with the isodensity contours generated by KDE method. Lines A and B 
represent the "axis" of the WISE strip occupied by BZBs and the directions of maximum 
gradient of the ratio between probability density of the BZBs to the WISE thermal sources. 
The inset histogram shows the normalized distribution of distances of BZBs sources (red) 
and WISE thermal sources (black) calculated along the line A (see Section [3]). Labels on the 
contours indicate the percentage of BL Lacs contained in each contour. 
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4. VER J0648+152: an unidentified TeV source 



The blazar position over tlie WBS can be used as diagnostic tool to identify extragalactic 
7-ray sources. We applied this to t he case of the TeV un identified source VER J0648+152 
recently discovered by VERITAS (lOng fc Paneaue 20101). t hat has been hkely associated 
with a blazar at low galactic latitudes ( IStephen et al. 2010l ). The blazar identific ation has 
been suggested on the basis of both the X-ray counterpart (IStephen et al. 20101 ) and the 
Fermi de tection of the 7-ra y source 1FGLJ0648. 8+1516 within the positional error circle of 
0.03 deg JAbdo et al. 201oh . 



We considered the positional circle of radius 108" = 0.03 deg, corresponding to the 
error on the position of the Fermi source and we report the infrared colors for all the WISE 
sources over the [3.4]- [4.6]- [12] /zm color-color diagram. We found that the closest WISE 
object (20") to the Fermi position lies out of the WBS in the area of the [3.4]-[4.6]-[12] /xm 
color-color mostly dominated by stars and elliptical galaxies. However, the second closest 
WISE source (29") lie in the region of the WBS where the TBLs are mostly concentrated 
(see Figure [3]), correspondin g to the position of t he X-ray counterpart within the error box 
of XMM-Newton (i.e., 7". 2) ( IStephen et al. 20101 ). In Figure[3l the open squares indicate all 
the sources associated with the VER J0648+152 unidentified object within the error box of 
the Fermi position, while the filled square points to corresponding the WISE source. The 
association of this WISE source strongly supports its blazar classification and 
given its position on the WBS, we suggest it is a new TBL. 



5. Summary and Discussion 

The most striking result of our investigation is that in the [3.4]- [4.6]- [12] /im color-color 
diagram blazars lie in a distinct region (WBS) with respect to the other extragalactic sources. 
This is due to the fact that blazars are dominated by non-thermal emission, while all the 
other extragalactic classes of galaxies by thermal radiation at infrared frequencies. This is 
strongly supported by the infrared colors of the WBS that are in agreement with a power-law 
model for the IR spectrum of the BZBs (see also Figure [2]). 

A similar att empt has been previously performed using the 2MASS archival data (e.g., 
Chen et al. 20051 ) , to investigate the behavior of blazars with respect to the normal galaxies 
in the J-H-K color-color diagram. However, our new approach has two main advantages: 
first, because we are using mid-IR wavelengths where the stellar photospheric contributions 
do not dominate galaxy colors, the blazar population covers a distinct region in the [3.4]- [4.6]- 
[12] /im color-color plot, that separates them from the other extragalactic sources; second. 
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because the Galactic absorption at the WISE wavelengths is negligible. 

The discovery of the WBS has several implications. When the WISE data over the 
whole sky will be available it will be possible place all sources over the [3.4]-[4.6]-[12] fim 
color-color diagram and those unidentified that lie in the WBS could be considered blazar 
candidates. Multifrequency observations can be used to rule out this hypothesis; this will 
be also simplified given the good astrometric precision of WISE data. 

An additional use of the regions in the WBS is that they provide diagnostic tool for the 
blazar associations. As previously discussed in Section |2l in our analysis we excluded the 
multiple associations between a WISE source and a blazar in the ROMA-BZCAT. However, 
our investigation shows that if at the position of the ROMA-BZCAT source within few 
arcseconds there is a correspondent source in the WISE blazar Strip, this could be identified 
and associated with the blazar. 

A samilar procedure could be applied to solve the problem of the unidentified 7-ray 
sources (see Section H]). In the 7-ray sky most of the extragalactic 7-ray sources are blazars, 
in particular those that are emitting at GeV-TeV energies, where the HBLs dominate. We 
showed the use of the WBS as diagnos tic tool for the case o f VER J0648+152 recently 
detected at TeV energie s by VERITAS 



Ong &: Paneque 2010l ). and already present in the 



1st year Fermi catalog ( Abdo et al. 2010[ ). Our analysis of the infrared colors for all the 



WISE sources in the Fermi positional error circle strongly support its blazar classification 



fIStephen et al. 20101 ). Its position on the WBS suggest that VER J0648+152 is a new TBI. 



Finally, we note that the region of the WBS covered by the TBLs is well defined. 
Then, when the whole WISE data is released this region could be used to identify new 
TeV candidates for future observations. A more accurate analysis on the possible use and 
developments of the WBS as a diagnostic tool to classify blazars as w ell as the unidentified 7 



ray sources will be pr esented in more details in a forthcoming paper ( iD'Abrusco et al. 2011 



Massaro et al. 2011bl ) 
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